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Abstract 

We study analytically and numerically the decay of a metastable phase 
in (2+l)-dimensional classical scalar field theory coupled to a heat bath, 
which is equivalent to two-dimensional Euclidean quantum field theory at 
zero temperature. By a numerical simulation we obtain the nucleation 
barrier as a function of the parameters of the potential, and compare it 
to the theoretical prediction from the bounce (critical bubble) calculation. 
We find the nucleation barrier to be accurately predicted by theory using 
the bounce configuration obtained from the tree-level ("classical") effective 
action. Within the range of parameters probed, we found that using the 
bounce derived from the one-loop effective action requires an unnaturally 
large prefactor to match the lattice results. Deviations from the tree-level 
prediction are seen in the regime where loop corrections would be expected 
to become important. 



1 Introduction 



In the early eighteenth century, Gabriel Daniel Fahrenheit [|IJ noticed that pure water 
could be cooled well below 32° on his newly invented temperature scale, and still remain 
a liquid. It was not until 1935, however, that Becker and Doring || gave a quantitative 
nucleation-theoretic treatment of droplet formation for fluid systems. Phenomenological 
field theory treatments were developed by Cahn and Hilliard and by Langer || 
within the context of a coarse-grained Ginzburg-Landau model. In quantum field theory, 
the study of metastable vacuum decay was initiated by Voloshin, Kobzarev, and Okun 
[H], and put onto firm theoretical ground by Coleman and Callan ||. The realm of 
applicability of homogeneous nucleation theory is extremely wide, from vapors, liquids, 
and solutions to metals, polymers, and glasses [0]. Recently, the possibility that first- 
order phase transitions occurred in the early Universe has generated a great deal of 
interest in the metastable decay of the vacuum itself. Well known examples are inflation 
[§], the electroweak phase transition ||, and the quark-hadron phase transition flO| . 



Surprisingly, however, classical nucleation theory has yet to receive clear experimental 
verification |TT[] : most experimental systems have complicated features that are ignored 
when formulating a Ginzburg-Landau description to which the classical theory could 
be applied. Furthermore, even when a reliable field-theoretic description is known, it 
is unclear how to accurately obtain the coarse-grained free-energy functional that is 
required for the critical-bubble calculation ]rj| . In cosmological applications, it has been 
customary to use the one-loop effective potential in the calculation of the nucleation 
barrier ||, |L3| . However, this procedure has been questioned, since Vi_i oop only reflects 



static properties of the theory, and hence may not be a good guide to the dynamics of 
out-of-equilibrium processes |14|- Also, the incorporation of temperature corrections to 



the nucleation barrier has not been carefully considered [15[ . 

Given the universality of the topic it would seem desirable to test the classical 
theory within its area of validity by performing a numerical simulation of some simple 
model, for which we can measure the effective nucleation barrier. This barrier would 



then be compared with theoretical predictions. In recent work |I7| we carried out this 
programme for a one- dimensional Ginzburg-Landau model, and showed that classical 
nucleation theory using the microscopic (lattice) Hamiltonian agrees well with numerical 
simulations. In this letter we will analyze the more realistic two-dimensional case, 
for which the lattice Hamiltonian is cutoff-dependent, and hence it cannot be used to 
determine the nucleation barrier. 

The system we study is classical thermal (2+l)-dimensional scalar field theory, which 
is closely related to zero-temperature quantum scalar field theory in two-dimensional 
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Euclidean spacetime. Replacing the temperature T with the quantum of action h, the 
equilibrium correlation functions of the classical theory map exactly onto the Green's 
functions of the quantum theory, and the free energy maps onto the effective action. 
Since we will find it convenient to mix the classical and quantum terminologies, we give a 
dictionary in Table 1. Even though the classical nucleation rate in Langevin time seems 
like a non-equilibrium quantity with no quantum analogue, classical nucleation theory 
identifies it (up to a dynamical factor) with the imaginary part of the equilibrium free 
energy density of the metastable state [[|. Under T — ^ h this maps onto the quantum 
tunnelling rate, which is given in terms of the imaginary part of the energy density of 
the false vacuum |J . Therefore, one could also think of our simulation as an independent 
way of calculating the barrier to zero temperature quantum tunnelling. In the classical 
language, we will show that the barrier obtained from the tree-level free energy (the 
Hamiltonian), ignoring entropic (loop) corrections, is in excellent agreement with the 
numerical results. In the quantum language, what we will show is that the tree-level 
(classical) action yields the correct bounce action in the Callan and Coleman calculation. 
Naturally this result may be modified when the scalar field interacts with other fields 

There has been only one previous numerical study of nucleation, that of Vails and 



Mazenko | I6|| , which found that the theory gave a very poor estimate of the nucleation 
barrier in a two-dimensional Ginzburg-Landau model coupled to a heat bath. There are 
several reasons |I7[ why our results should differ from theirs, the dominant one being 
that two-dimensional classical thermal field theory suffers from ultra-violet divergences 
which produce severe lattice-spacing dependence of many measured quantities, including 
nucleation rates. Continuum results can be obtained, as we will show, by introducing 
cutoff-dependent counterterms and renormalizing. Since Vails and Mazenko did not 
explicitly introduce counterterms, it is possible that the continuum theory they were 
simulating was different from the one for which they were calculating theoretical decay 
rates. 

We will now review the "critical bubble" or "bounce" method of calculating the decay 
rate, before proceeding to test it numerically. 



2 The Theoretical Results 

Consider a classical (2+l)-dimensional scalar field theory with an asymmetric double- well 
potential, at finite temperature. By virtue of its relation to a two-dimensional quantum 
field theory, mentioned above, we expect this system to be formally ill-defined because 



2 



Classical Statistical Mechanics 


Euclidean Quantum Field Theory 


Microscopic (Lattice) Hamiltonian 

-niatt = J rtlatt « X 


Bare Action 

*Jbare J *-bare " 


Lattice Energy Density 

H-latt = 2^1att| V0r + ^latt(0; «) 


Bare Lagrangian Density 

f 1 V IV7J.I2 i \r I A „\ 
Mwe — ^barel V<p| + Vbarel'P, a ) 


Free Energy 
F = /^d 2 a; 


Effective Action 

S C S = I Ceft d 2 X 


Free Energy Density 
F=\Z F \V<P\ 2 + V F (<P) 


Effective Lagrangian Density 
£ eff = iZ cff |V0| 2 + Kfr(0) 


oo 

F(0,T) = H{cj>) + ]TAF n _ loop (0)T n 


CO 

Ses(<f>, h) = Sdass(0) + H AS^loop^) ft™ 


Entropic corrections 


action Quantum corrections 


, " V 

OO 

V F ((f>,T) = U((f>) + AK- loop (0)T" 

?1=1 


, A , 

CO 

V eS (<P,h) = U((j>) + £ AK-ioo P (0)^ n 

Classical 
potential 



Table 1: Dictionary of equivalent terms for classical thermal and quantum mechanical 
two-dimensional field theory on a lattice with spacing a. 

of ultra-violet divergences in the correlation functions. For the moment we therefore 
imagine integrating out the high-momentum physics, and assume that the low-momentum 
dynamics is described by some coarse-grained free-energy functional, 

^F dyn [0] = ly {±Z dyn |V0| 2 + V Ayn {cf>)}d 2 x . (2.1) 

We assume that Vd yn (</>) has a local minimum at cf) = 0, separated by a barrier from a 
lower global minimum at some positive (perhaps infinite) <ft. 

In general one does not know what F& yn should be. Langer [0] states that it should 
be a free energy coarse-grained up to the correlation length (but no further). We will 
investigate two candidates for -Fd yn , the tree-level free energy or Hamiltonian H (the 
classical action, in the quantum language) and the 1-loop free energy (1-loop effective 
action) (see Table 1). 

We will study nucleation for potentials of the general form 

U(<j>) = \m 2 4> 2 - \g<f + ±h<p A . (2.2) 
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This seems to give a 4-parameter family of theories, but in fact by rescaling the field and 
distance, we can reduce it to a functional involving only two (dimensionless) parameters, 
A and 9: 

h m = \\ + & " ^ + ^ A04 i ^ ' (2 ' 3) 

where 

„ g 2 T , , hm 2 . . 

9 = y and A = — — 2.4 
m 4 g 2 v ' 

The dimensionless parameter 9 plays the role of temperature, while A parameterizes the 

different potentials. The potential has one minimum at <fi = 0, and for < A < | there 

is a second, lower minimum at positive 0, separated from the metastable minimum by a 

maximum at max = (3/2 A) (l — ^1 — 8A/3^. 

For such a system, we expect Fd yn to be a functional parameterized by A and 9, such 

that if the system starts out with thermal expectation (<fi(x)) = at finite temperature 

9 then it will make a thermally activated transition to the global minimum with a 

characteristic rate per unit area 

r = Aexp(-B/9) , (2.5) 

where A, the prefactor, is a function of A, 9, and the coupling to the heat bath rj (see 
Sect. 3). B is the nucleation barrier, which, according to the classical nucleation or bounce 
calculation [|, |[, is the energy of the "critical bubble" or "bounce" saddle point of Fd yn . 
The critical bubble configuration (f>B (x) is given by the static (circularly symmetric) 
solution to the equation of motion, 

^r + -^r = ^ yn ^(<M , (2.6) 

with the boundary conditions that the field be regular at the origin (0'#(O) = 0) and in 
the false vacuum at spatial infinity {(f) b (oo) =0). If one thinks of r as being "time", this 
equation describes the motion of a particle with position (f> B in the potential — Vd yn (0), 
with a frictional resistance that is inversely proportional to the time. In general there is 
no closed- form solution to Eq. (|2.6|) , but it can be solved numerically by guessing a value 
for 0b(O) and evolving forward in r. If the solution overshoots and 4>b becomes negative 
then the initial value was too high; if dips/dr becomes positive then it was too low. In 
this way we converge on the correct solution, and its energy is the nucleation barrier B. 

Since in general F& yn may depend on 9, we may write B = B + B x 9 + B 2 9 2 H , and 

so Eq. ( |2.5| ) becomes 

r = A exp(-Bi) exp [- (B /9 + B 2 9 + • • ■)] . (2.7) 
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As was mentioned above, the tree-level candidate for Fd jn is the Hamiltonian H. 
This is independent of 9, so B will just be Bo, and we can calculate it by the method 
described above, setting Vd yn = U. Another possibility is that F dyn is H plus some 
entropic corrections. For example, in studies of cosmological phase transitions it is 
generally assumed to be the real part of the 1-loop free energy. In this case Bi, B 2 etc 
are non-zero, but as 9 — > the entropic corrections disappear, so B is still given by the 
tree-level (-Fdyn = H) barrier prediction. The lowest order correction B\ only modifies the 
prefactor, not the measured barrier. Thus entropic corrections will primarily be visible 
through the effects of the B 2 term, since we have little information about the theoretical 
value of A beyond the fact that we expect it to be of order 1 on dimensional grounds. 
It will turn out that we can confirm that Bq is given by the tree-level bounce, but our 
observations cover too narrow a range of 9 for the B 2 predicted by the 1-loop free energy 
to be visible, so we are unable to say definitively whether it gives the right entropic 
corrections. 

We expect the barrier to go to infinity when the two minima are degenerate, at A = |, 
and to go to some finite value as A —>■ 0. In fact, A can take on negative values and 
everything still works exactly as described above. There is no longer a global minimum 
to the potential, but there is still a metastable state which decays by surmounting a 
barrier that can be calculated by finding the critical bubble. We will restrict ourselves 
to studying A in the range to 0.32, since negative A is not physically relevant, and 
as A — > | the barrier becomes large, requiring high temperatures to induce nucleation, 
which invalidates the loop expansion we use in the next section. 



3 The Lattice Formulation 

In order to test the critical bubble theory, we need to simulate the classical dynamics of 
the (2+l)-dimensional scalar field theory in contact with a heat bath at temperature 9. 
This may be done by evolution of the stochastic Langevin equation, 

d 2 (p 2 <90 dV l&tt . . 

^T = ^ t V0-,-- — + £(M) , (3.1) 

where rj is the viscosity coefficient, and £ is the stochastic noise with vanishing mean, 
related to r\ by the fluctuation-dissipation theorem, 

(£(£, t)f (f , t')) = 2 V 95{t - t')5 2 {x - x' ) . (3.2) 



In writing Eqns. (|3.1| ) and (|3.2| ), we implicitly assumed that the system is Markovian, i.e. 



the correlation time scale for the noise is much smaller than the typical relaxation time 
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for the system, which is the inverse of the oscillation frequency around the metastable 
equilibrium. 

In principle it now seems straightforward to put the field theory on a two-dimensional 
lattice, and solve Eq. numerically. However, ( |3.1|) is expressed in terms of the lattice 
(bare) parameters, so in order to make contact with the calculation of Sect. 2 we must 
explicitly construct a lattice potential that has a definite continuum limit (lattice spacing 
a — > 0), and have some way of specifying which continuum theory it corresponds to. Since 
the theory has exactly the same divergences as a two-dimensional Euclidean quantum field 
theory, the lattice potential will have to contain lattice-spacing-dependent counterterms 
in order to give a good continuum limit. (For a review, see Parisi jn|, Ch. 5). 

Our strategy will be to use continuum field theory with a hard momentum cutoff 
A ~ a -1 as our guide. We can easily calculate the free energy (effective action) for such a 
theory with a given bare potential by a loop expansion, which is an expansion in powers 
of 9 (h in the quantum theory). This will tell us what counterterms are needed to cancel 
the cutoff dependence of the lattice theory. We will only have to calculate to one loop, 
since there are no divergent graphs beyond one loop. 

We studied the lattice action defined by (see Table 1): 

Matt(0) = C/(0)+Vct(0) , 

W) = |0 2 -|0 3 + ^A0 4 , ( 3 - 3 ) 

V ct (<P) = --^ln(aM 1 (A,^)0+^ln(aM 2 (A,^))0 2 . 
Air Sir 



By the standard calculation |19| for a scalar quantum field theory in two Euclidean 
dimensions with classical potential U((f>) and ultraviolet cutoff A, we find that the 1-loop 
effective potential is 

M r\ 2 h / T7"\ 

Performing the integral and dropping terms independent of 0, we find for large A, 

V F = U + V ct + ^-(U"\n(A 2 ) + U"-U"HU")) . (3.5) 

The counterterms cancel the A-dependence, so we expect the theory with lattice potential 
( |3.3|) to have a good continuum limit. What effective potential will this continuum theory 
have? Eq. ( |3.5|) was calculated for a theory with momentum cutoff A to simulate the 
effects of the lattice. It does not tell us what finite parts will be left from the cancellation 
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of the cutoff-dependence in the actual lattice theory. Thus we can only say that we 
expect the lattice theory defined by (|3.3| ) to have one-loop effective potential 

Vf = U(<f>) + ^ (U"{(j>) - U"{4>) HU"{4>))) + f{M x )4> + <?(M 2 )0 2 . (3.6) 

Mi and M 2 are to be fixed by a renormalization condition, but since the functions / 
and g are unknown, this will have to be done numerically for the lattice theory. Our 
renormalization condition was 

f(M 1 ) = g(M 3 ) = . (3.7) 

It only remains to calculate the 1-loop effective field normalization Z (see Table 1). 
It is given in terms of the one-particle- irreducible Green function T^(p), and is cutoff- 
independent in two dimensions. Evaluating the relevant Feynman diagram, 

dT ( V 9 

To summarize, the lattice theory ( |3.3| ) with renormalization conditions ( |3.7| ) will have 
one-loop free energy density (effective action) 



m 2 + u(<p) + ^-(u"(4,) - u'\4>) Hu"m) (3-9) 



^=( 1 -iL)i v ^ + f/ w + ^v 

U" is negative between the points of inflection of U, so T is complex. Weinberg and Wu 



21j| have suggested that the real part of this is physically meaningful, so when using T 
in bounce calculations, we take the real part. 

All these calculations are valid to one loop. Without explicitly evaluating the higher 
corrections, we can estimate that they will be valid for 9 <C 8n, at least near the minima 
of the potential. There are no additional divergences at two loops, but finite corrections 
to J 7 , proportional to 8 2 , will appear. This puts a limit on how close we can push A to 
the value 1/3, since the barrier becomes large as the the two minima in U(4>) become 
degenerate, requiring high temperatures in order to observe nucleations on the lattice. 



4 The Numerical Analysis and Results 

We now have a way to test the theoretical calculation of Sect. 2. For each value of A we 
measure the nucleation barrier on the lattice by measuring nucleation rates for a range 
of values of temperature 9. For each A and 9 this involves two stages. Firstly, we ensure 
that we are looking at the right continuum theory, by running simulations with lattice 



7 



potential (|3.3| ), and varying Mi and M 2 until ( |3.7p is obeyed. We then know that we 
have a theory whose free energy is given by (|3.9|) , to first order in 6/8tt. The second stage 
involves measuring the nucleation rate for this theory. In both stages we ensure that we 
are studying the continuum limit by reducing the lattice spacing a until the results become 
independent of a. We checked the dependence of our results on the lattice length, the 
time step, the random number generator and the random number seed. Within the limits 
of our numerical accuracy we found that the lattice approximation correctly describes the 
continuum field theory. We chose 77 = 1 as the viscosity in all simulations. In principle, 
knowledge of the physical nature of the heat bath would enable one to calculate the 
value of 1], as well as additional nonlinear or nonlocal dissipative terms in the Langevin 
equation. 

The easiest way to impose the renormalization conditions was to choose Mi and M 2 
such that the measured true and false vacua on the lattice coincided with the minima 
of the free energy (|3.9|) . This required running a simulation of Q3.1|) on a relatively 
small lattice (L = 10), with initial conditions that caused the field to settle down in 
the appropriate vacuum state. It was then a simple matter to measure the thermal 
expectation value of the field. We found that the renormalization condition /(Mi) = 
(i.e. ((p) = in the metastable state) was obeyed for Mi = l/(2n), in the sense that 
with this choice the measured metastable average field was always much less than the 
peak max . The condition g(M 2 ) = was obeyed for M 2 = l/(27r), in the sense that 
the measured true vacuum was within a few percent of the value predicted by ( |3.9| ). 
Obviously this method could not fix M2 in the A = theory, since the true vacuum is 
then at infinite 0. However in that case there is no cutoff-dependent 2 term in the free 
energy, and g(M 2 ) = 0. 

To measure the nucleation rates, we evolved the Langevin equation on square lattices 
of size L = 20, 40, 60 or 120, with lattice spacing a = 0.5, using a leapfrog algorithm with 
time step of 0.05. For each A and 9 we performed several hundred simulations. In each 



one we started the system in the metastable state, evolved the Langevin equation ( pTT 



forward in time, measuring the time t that elapsed before it escaped from the metastable 
region by nucleating a growing bubble of the stable state. From this data we constructed 
a frequency histogram for nucleation times from which T(X,6), the nucleation rate per 
unit area, could be read off by fitting to the expected form: Prob(t) oc exp(— tTL 2 ) for 
a lattice of size L. At low temperatures T gets small ( [2.5|) , but we could compensate for 
this by increasing L. Thus with a range of lattice sizes we were able to observe a wider 
range of temperatures than would have been possible using only one lattice. We were still 
constrained to a fairly narrow temperature range, since the simulation ran much more 
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slowly on larger lattices. 

In order to determine when nucleation had occurred, we calculated a smoothed field 
by averaging over square blocks of size AL. The system was considered to have escaped 
from metastability when one or more of the blocks achieved an average field value greater 
than the peak max (see Sect. 2). For each value of A we chose the smallest AL such 
that when one block was converted, the rest of the lattice would always follow. For 
A = 0, 0.1, 0.2 we set AL = 5; For A = 0.25, 0.3, 0.31 we set AL = 10; for A = 0.32 we set 
AL = 20. 

Our results are summarized in Fig. 1, 2, and 3. 

In Fig. 1 we show a logarithmic plot of the inverse nucleation rate per unit area r _1 , 
as a function of inverse temperature, for different values of the potential parameter A. 
From Eq. for the decay rate we expect that ln(r -1 ) = \n(A- 1 ) + B l + B (X)/9 + B 2 9. 
The excellent straight line fits for A < 0.3 indicate that the classical theory that gave 
rise to Eq. (|2.7p is qualitatively correct in this regime. There is no sign of a B 2 9 term. 
We conclude that the nucleation barrier is well approximated by the form B = B Q + B\9 
for the ranges of temperature probed for A < 0.3. For A = 0.31, 0.32 there are signs of 
a B 2 term, so entropic corrections to the barrier are becoming important, as one would 
expect at high 9. 

In Fig. 2 we confront the numerical barrier measurements from Fig. 1 with the 
theoretical prediction for the nucleation barrier B given in Sect. 2, which involves taking 
Vd yn (0) to be U (</)), the potential appearing in the tree- level free energy (Hamiltonian, or 
classical action). For A < 0.3 we obtain almost perfect quantitative agreement between 
theory and numerical experiment, confirming that B is accurately given by the action 
of the tree-level bounce. As one might expect from the weakness of their straight line fits 
in Fig. 1, the barriers for A = 0.31, 0.32 are not well predicted by the tree-level bounce. 

In Fig. 3 we plot the values of the prefactor A that would fit the data of Fig. 1, using 
the tree-level bounce and the 1-loop bounce. For the tree-level points, we fitted Eq. ( |2.7|) 
to the data of Fig. 1 with Bi = 0, and plotted the value of A' 1 for each value of A. 
For the 1-loop points we calculated B(\, 9) using the method of Sect. 2 with Fd yn set 
to the real part of the 1-loop free energy ( |3.9|) . (This was only possible for 0^10. For 
larger 9 no solution could be found, so 1-loop barrier predictions could not be given for 
A = 0.31, 0.32, whose rate measurements were all at larger 9.) By fitting to the form 
B = £>o(A) + Bi(X)9 we obtained B (which agreed with the tree-level prediction, as 
expected) and B\. For each A we then found A by fitting Eq. ( |2.7| ) to the data of Fig. 1. 

Much larger prefactors were required to fit the 1-loop predictions to the data than 
for the tree- level predictions. Since A is expected to be of order 1, this provides 
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circumstantial evidence that the one-loop bounce is not the relevant quantity. For a 
definitive answer, a more powerful computer would be needed, allowing a larger range of 
6 to be probed for each A, so that the value of B 2 could be measured. 

The fact that the results fit the theory less well as A approaches | should not be 
surprising. Since the barrier is diverging in this region, we are forced to use high 9 
in order to see any nucleations, and so we expect entropic corrections to the barrier 
to become important. The deviations begin at A = 0.31, where temperatures in the 
range 12 to 20 are investigated. This is approaching the regime 9^8n, where the 
loop approximation breaks down, so we expect significant ^-dependent corrections to 
the theoretical nucleation barrier. In other words, we expect B 2 and higher terms to 
become large. It is perhaps more surprising that we do not see deviations for A = 0.3, 
where a similar temperature range was explored. 



5 Conclusions 

We have shown that, for a single scalar field in two space dimensions, classical nucleation 
theory correctly predicts the nucleation barrier to within a few percent, when the 
Hamiltonian, i.e. the tree-level free energy, is used in the calculation. Equivalently (see 
Table 1), the quantum bounce calculation correctly predicts the barrier, as long as the 
tree-level ("classical") action is used. 

Obviously this statement only makes sense for the parameter range where a loop 
expansion is possible, and, as expected on the basis of the calculations of Sect. 3, we 
found that it broke down for large loop parameter, 6/8^^1. The interesting thing is 
that the tree-level action gave correct predictions for 8/8tc as large as 0.5, a regime where 
one might have expected significant loop corrections. 

We were unable to make any definitive statement about the 1-loop bounce (the bounce 
solution obtained by using the 1-loop effective potential). Fig. 3 shows that it requires 
an unnaturally large prefactor, but the temperature range we probed was too narrow to 
see if its predicted higher-order corrections were present. With the aid of more powerful 
computing resources it should be possible to settle this question. 

There are many other directions in which this work can be extended. E. Weinberg 



2^] has recently argued that a truncated 1-loop effective potential correctly describes 
nucleation in certain radiative symmetry breaking models, and it would be interesting 
to study this numerically. Another issue is finding the dependence of the prefactor not 
only on the temperature but also on the coupling to the heat bath. We have previously 
studied the dependence of the thermal kink-antikink pair creation on the coupling to 
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the bath [20], finding good agreement with Kramer's celebrated result |[23|| . Perhaps of 



more relevance to current work on out-of-equilibrium processes in the early Universe, we 
could study the dynamics of weakly first-order transitions. Recent work has claimed that 
for sufficiently weak transitions the usual bubble nucleation mechanism discussed here 
must be modified. Instead, it has been suggested that the transition may evolve through 
the nucleation and subsequent percolation of sub-critical bubbles in such a way as to 
resemble an emulsion of phases as the system is cooled below the critical temperature 
?M. We are presently investigating these questions. 
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Figure Captions 



Figure 1: Numerically measured inverse nucleation rate T as a function of inverse 
temperature, with straight line fits, for a range of values of the dimensionless 4 coupling 
A. From the right: A = 0, 0.1, 0.2, 0.25, 0.3, 0.31, 0.32. 

Figure 2: Nucleation barrier as a function of A. The points are the slopes of the lines 
in Fig. 1. The line is the critical bubble calculation of Sect. 2 at tree level, i.e. with 
Z dyn = 1 and V dyIl = U. 

Figure 3: Rate prefactor A^ 1 as a function of A. Bottom set of points is for the tree-level 
prediction (see Eq. (|2.5|) ) fitted to the data of Fig. 1. Top set of points is for the 1-loop 
prediction (see Eq. (|2.7|) ). We excluded A = 0.31, 0.32 because their error bars were too 
large. 
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